arXiv:1505.07473vl [hep-ph] 27 May 2015 


A 125 GeV Fat Higgs at large tan 13 


Arjun Menon®’^ and Nirmal Raj^ 

Institute of Mathematical Sciences, Taramani, Chennai, 600113, India 
’^Institute of Theoretical Science, University of Oregon, Eugene, OR 97405 

May 29, 2015 


Abstract 

In this article we study the viability of regions of large tan /3 within the framework 
of Fat Higgs/A-SUSY Models. We compute the one-loop effective potential to find the 
corrections to the Higgs boson mass due to the heavy non-standard Higgs bosons. As the 
tree level contribution to the Higgs boson mass is suppressed at large tan/I, these one- 
loop corrections are crucial to raising the Higgs boson mass to the measured LHC value. 
By raising the Higgsino and singlino mass parameters, typical electroweak precision 
constraints can also be avoided. We illustrate these new regions of Fat Higgs/A-SUSY 
parameter space by finding regions of large tan S that are consistent with all experimental 
constraints including direct dark matter detection experiments, relic density limits and 
the invisible decay width of the Z boson. We find that there exist regions around 
A = 1.25, tan/? = 50 and a uniform psuedo-scalar 4 TeV ^ Ma ~ 8 TeV which are 
consistent will all present phenomenological constraints. In this region the dark matter 
relic abundance and direct detection limits are satisfied by a lightest neutralino that is 
mostly bino or singlino. As an interesting aside we also find a region of low tan /3 and 
small singlino mass parameter where a well-tempered neutralino avoids all cosmological 
and direct detection constraints. 



1 Introduction 


Weak scale supersymmetry (SUSY) remains a popular and elegant solution to the hierar¬ 
chy problem of the Standard Model (SM) [1]. It provides a natural means to stabilize the 
electroweak scale against large quadratic corrections from higher scales. The fact that 
the Higgs boson was discovered to be light mm, in accordance with SUSY’s predictions, 
encourages us to continue our search for signals of SUSY at the TeV scale. 

In the Minimal Supersymmetric extension of the Standard Model (MSSM), the tree- 
level Higgs quartic couplings are fixed to be the gauge couplings which leads to the 
tree-level Higgs boson mass that is below that of the Z boson. Therefore, raising the 
Higgs boson mass to the observed value of 125 GeV at the LHC mm requires large 
corrections due to a heavy stop sector m 0 0 0 El El nni dn Ea da HU . However, heavy 
stops lead to large corrections to the up-type soft SUSY breaking Higgs squared mass 
parameter which in turn leads to a large correction to the electroweak symmetry breaking 
(EWSB) condition. A delicate cancellation between these corrections and the Higgsino 
mass parameter is needed to stabilize the electroweak scale, which is generally considered 
unnatural. Therefore in the MSSM there exists a tension between the observed Higgs 
mass and the requirement that the model is natural. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) is the simplest 
extension of the MSSM that can address this tension. In the NMSSM, the Higgs sector 
is enlarged to include an extra gauge singlet that couples to the remaining MSSM Higgs 
doublets through a Yukawa coupling A. A contributes to the Higgs quartic at tree- 
level, and for large enough values, can raise the Higgs mass to the observed 125 GeV. 
Therefore the stops need not be too heavy, thereby improving the naturalness of the 
model. Moreover, in the general NMSSM (GNMSSM), an additional tadpole term for 
the gauge singlet can also facilitate EWSB m- 

Eor A > 0.7 at the weak scale, renormalization group (RG) evolution usually leads to 
this coupling developing a Landau pole below the GUT scale. Refs. US da da da EQi EH 
have provided explicit UV-completions for such low scale models, which we collectively 
call Eat Higgs models. Refs. [SaEaEllEaEaEaEaiZaBaEaEaEaEl! have 
studied the phenomenological implications of models with such large A couplings, which 
we collectively call A-SUSY models. For these models, they have found that the Higgs 
mass can easily be raised to the observed value while still keeping the spectrum natural. 
These studies have focused on a region of low tan/3 (= Vujvd, where Vu and Vd are the 
vacuum expectation values (VEVs) of the corresponding Higgs doublets) and large A 
because these regions were the most natural. 

In this paper, we study the possibility of raising the Higgs mass to 125 GeV in Fat 
Higgs/A-SUSY models at large tan/3. As the A^-proportional tree-level contribution to 
the Higgs quartic is suppressed at large tan/3, the one-loop induced radiative correc¬ 
tions are crucial in raising the Higgs boson mass to its observed value. Similar to the 
stop-induced corrections that are proportional to yf log{m‘~/Q'^) (where yt is the top 
Yukawa and is the stop mass scale), in Fat Higgs/A-SUSY models the dominant one- 
loop corrections are proportional to A^ log(M^/(5^) (where Ma is the scale of the non¬ 
standard Higgs bosons). Therefore these corrections are only relevant when A>1 and 
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the non-standard Higgs bosons are much heavier than the electroweak scale. The effect 
of radiative corrections in the NMSSM Higgs sector have been considered before |29l l35] . 
Ref. [29] focused on the most natural regions in the Scale-Invariant NMSSM, where it 
was found that these radiative corrections made a negligible contribution. In contrast, 
we show that at large tan /3, the A induced radiative corrections can significantly modify 
the allowed regions of parameter space. Unlike Refs. [221 EH EO], we also emphasize 
that electroweak precision constraints do not put a limit on tan /3. We point out that 
raising the Higgsino mass parameter //eff significantly weakens the electroweak precision 
constraints because the Higgsino component in the lightest neutralino is suppressed. The 
price of raising /Tefr is a slight increase in the tuning of the EWSB condition. To illustrate 
these effects in regions of large tan /3 we also impose constraints from Higgs decay prop¬ 
erties, direct dark matter detection experiments, the observed dark matter relic density 
and the invisible width of the Z boson. In particular, we find that direct dark matter 
detection experiments place strong limits on many regions of parameter space due to the 
large A coupling. We also show that these allowed pockets of parameter space are within 
the reach of the XENON IT experiment jnSj. 

This paper is organized as follows. In Sec. [2| we set up the theoretical aspects 
required for the phenomenology of our model. To motivate the sizes of various terms 
in the Fat Higgs/A-SUSY model, we present a “toy” high scale model where the fields 
have canonical mass dimensions in the electric theory. In addition, in this section we 
also compute the corrections to the Higgs quartic using the one-loop effective potential 
formalism, and discuss the Higgsino contributions to electroweak precision constraints 
and naturalness in the large tan/3 regime of the Fat Higgs/A-SUSY model. In Sec. 
we illustrate the impact of the formalism in Sec. by finding phenomenologically viable 
scenarios that can be probed at future experiments. In Sec. we conclude. 


2 Theoretical Setup 

In this section we first motivate the form taken by our superpotential by a simple dis¬ 
cussion of the sizes of various terms that can arise in Fat Higgs/A-SUSY models. In 
this discussion we assume that any exotic fields are much heavier than the electroweak 
scale. For the superpotential thus obtained, we present the Higgs potential at the tree 
level and analytically compute the one-loop corrections to the mass of the SM-like Higgs 
boson due to heavy non-standard Higgs fields, with special attention to the limit of large 
A and tan /?. In addition we discuss the naturalness of the large tan /3 regions of the Fat 
Higgs/A-SUSY models. We then discuss the reduced couplings of the SM-like Higgs to 
SM particles, which are constrained by LHC measurements of signal strengths. We end 
the section with a brief discussion of the neutralino sector with particular attention to 
electroweak precision observables. 
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2.1 Realizing low scale NMSSM with large A 

The GNMSSM with a large A at the weak scale implies that some of the Higgs fields 
are composite states. For example, in the minimal Fat Higgs scenario of Ref. [TB], all 
of the Higgs sector fields are composite, while in Refs. [m El] the MSSM Higgs fields 
are fundamental. For simplicity we will assume that at scales ^ 10 TeV, the only fields 
present in the Higgs sector are the SU{2)l doublets Hu, H^ and the gauge singlet 50 
The most general superpotential with this particle content (assuming R-parity) has 
the form [JS] |T5| 

WgNMSSM = WYukawa + ^SHuHd + 3 ^*^^ + l^^uHd + (1) 

where A, k are dimensionless coupling strengths; //' are supersymmetric mass terms; 

is a supersymmetric tadpole term of mass dimension 2 , and VVyukawa contains the 
standard MSSM Yukawa superpotential terms. The corresponding soft SUSY-breaking 
terms are 

- -^soft = -^Lft + + 

(^A\HuHdS + —kAkS^ + rn^HuHd + 2^5 h.c^ , ( 2 ) 

where corresponds to the standard MSSM soft SUSY-breaking terms. 
are the soft SUSY breaking Higgs squared mass terms and A\,Ai^ are the soft SUSY 
breaking trilinear terms. m|, m'^ are the CP-violating soft SUSY breaking squared mass 
terms and ^5 is the dimension-3 soft SUSY breaking term corresponding to 

A generic feature of most Fat Higgs/A-SUSY models is that the Yukawa coupling 
A > 0.7 at the TeV scale0 Due to its renormalization group (RG) evolution, A becomes 
stronger at higher scales and develops a Landau pole at the compositeness scale Ah, 
where Kp is assumed to be lower than the grand unification scale Mqut- In the deep 
IR, much below Ap, the magnetic theory of mesons (i.e. the Higgs superfields) is de¬ 
scribed by the interactions in Eqs. Q and Q. In the UV above Ap, some or all of 
the Higgs superfields are revealed to be composite states made up fundamental quarks 
whose interactions are described by some electric theory. 

If the quarks in the electric theory have the canonical mass dimension and all Higgs 
superfields are composite, then the k, //' terms in Eq. Q and their corresponding soft 
SUSY-breaking terms in Eq. Q are generated by marginal terms in the fundamental the¬ 
ory. This assumption is equivalent to saying that the confining dynamics only generates 
the A term while all other couplings arise from non-renormalizable interactions in the 
electric theory. For example, in the simplest Fat Higgs model [in]> the Higgs superfields 
in Eq. 0 are composite states of the quarks Ti in the electric theory. These quarks are 

^More exotic realizations typically require that additional Standard Model superfields are composite, in¬ 
troducing more exotic fields in the low energy theory [IllIllSTlEHI. 

^For Fat Higgs models that provide an existence proof of gauge coupling unification, see [igin]. 


3 



(3) 


charged under a confining SU{2)h gauge group, thereby leading to the identification 


5 ~ fgfb; 



(fits 

\f2n 




( Tif4 


The A term in Eq. 0 is dynamically generated by the Pfaffian of the mesons in the mag¬ 
netic theory. Naive dimensional analysis (NDA) |39 t IdO t BT] and canonical normalization 
of the fields ((TjTj) —>■ {ku /lead to the relations 


A(Aff) ~ dvr; 




I 

(47r)2Ao ^ ^ ’ 


k{Ah) ~ 


dvrAo 
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(4) 


where m and Aq are parameters in the electric superpotential given by 


bVelectric — mf^TQ + ^ (f^fo) + ^ 


An 


An - 


T 1 T 3 T2r4 - rir4 


+ 


A3 


iT5nf. 


(5) 


The couplings y,y',y", in the above equation, need not be 0(1) numbers because Aq is 
just a generic scale used to parameterize the mass dimension of each of these operators. 

Eq. 1^ gives us a definition of Ah- it is the scale at which the size of A is dvr. Using 
this definition, we can then estimate the size of the other parameters at the weak scale 
from their RG evolution. In determining Ah, we also account for the effects of the SM 
Yukawa couplings using the renormalization group equations (RGEs) in Ref. [33]. Having 
estimated the NMSSM parameters at the scale Ah using Eq. we run them down to 
the TeV scale by solving the RGEs and find that they decrease with decreasing scale. 
This running behavior has two important implications for our model: 

1. Eq. Q implies k{Ah) <C 0(1). Run down to a renormalization scale Q = O(TeV), 
we expect k to be quickly suppressed due to the contribution of A to its running. This 
suppression is illustrated in Fig. [^a), where we plot k at Q = Mz as a function of A at 
Q = Mz, setting tan/3 = 50. These curves were obtained by first running X(Q = Mz) 
up to determine Ah, then setting k(Q = Ah) to different values < 1, and finally running 
K down to Q = Mz- We checked that the running of A is insensitive to k for these sizes 
of K. The red curve corresponds to k{Ah) = 1 and the green curve to k{Ah) = 0.5. 
As expected from the RG running, smaller values of k{Ah) result in smaller values of 
k{Q = Mz). 

A larger A implies a Landau pole at a lower scale. Therefore, Ah is closer to the 
electroweak scale for larger values of A, which in turn weakens the suppression of k as 
it runs down from Ah to Mz- This is why k is an increasing function of A in Fig. 0a). 
From the plot, we infer that for k(Ah) < 1, the size of k at the weak scale is suppressed 
by at least an order of magnitude. The implication of this suppression is that we can 
consistently neglect the effects of n in our TeV-scale phenomenology. Therefore, for the 
rest of this paper we will take k = 0. 
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Figure 1: (a): k as a function of A at the scale Q = Mz, obtained by fixing k at the scale Ah 
and then evolving it down with RGBs. The red (green) curve corresponds to k{Ah) = 1(0.5). 
(b): /i' as a function of A at the scale Q = Mz, obtained in a manner analogous to (a). The red 
(green) curve corresponds to ^i'[Ah) = 1(0.5) TeV. In both plots we set tan/3 = 50. See text for 
details of their behavior. 


2. As compared to k, is only suppressed by an 0(1) number when it is run down 
from Q = Ah and Q = Mz- This difference between values of and k can be understood 
from their /3-function dependences. Using their one-loop /3-functions in Ref. [13] we find 

... 

k{Ah) VF'(Arf)y ■ ^ ^ 

We check this by determining /i'(Q = Mz) as a function of X{Q = Mz) in a manner 
analogous to the determination of k{Q = Mz) above. Our results are shown in Fig. [^b), 
where the red (green) curve corresponds to /i'(Ajy) = 1(0.5) TeV, with tan/3 = 50. We 
see that //'(Ajy) is suppressed at the electroweak scale by at most a factor of 5. Hence 
^'{Mz) can be of the size of the electoweak scale. Such a size results from reasonable 
values of the scale Aq. For instance, to obtain fi'{Q = Mz) ~ 1 TeV at tan/3 = 50, we 
computed from Eq. Hthat Aq < 10^® GeV for X{Q = Mz) > 0.7. Similarly, the Higgsino 
mass parameter n and the tadpole term can also be the size of the electroweak scale. 

We can now write down our low energy superpotential below the scale Ah'- 

^NMSSM = kVvukawa + XSHuHd + 2 ^^' (7) 

The associated soft-SUSY breaking potential is 

“ -^soft = + rn\\S\^ 

+ ^XA\SHuHd + rn^HuHd + T isS + h.c.^ . (8) 

We have redefined the singlet chiral superfield. S' —t' S' — /r, to remove the /r term in 
the superpotential. In general, the associated soft term cannot be absorbed into A\ 
simultaneously. Eqs. 0 and Eq. ([^ constitute all the parameters treated in the rest of 
this article. 
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2.2 Higgs Sector 

2.2.1 Tree level 

At the tree level, the Higgs potential is given by 

V^^^,, = Vf + Vd + Vs (9) 

where 


Vf 

Vd 

Vs 


|A - ilX) + ^^'S + CfI" + |A5|2 + |il,|2), 




+ ^ cos2 ew + HlH^V > (10) 




(AA,(il+il- - HlH^^)S + ml{H+H- - + h.c.), 


Hu = (-H^, -Hm), Hrf = (-Hrf) ^,7)’ = 51 + 52 is the weak mixing angle. After 

electroweak symmetry breaking we can expand the Higgs fields in terms of the CP-even 
fields {hu,h^,h^), the CP-odd fields the charged Higgs bosons and the 

Goldstone bosons (G^,G°); 



where v ~ 174 GeV is the VEV of EWSB, = sin/3, cp = cos j3 and s = (5). Expanding 
the potential about the minimum at Vi = {vu,Vd,s), we can find the tree-level tadpole 
terms 


Qytree 

octree _ ^ Higgs 

- dcP, 


f^i} 


( 12 ) 


where (pj = S). We can then solve for the soft squared masses 

by setting each pP®® = 0. Substituting these masses into the second order derivatives of 
the Higgs potential and neglecting CP-violating effects, we obtain the following tree-level 
CP-even Higgs mass matrix in the basis (/i^, h^, h^). 


(-l^lr)ii = ^ {^h)i 2 — (2A^u^ — Mz)spcs — r; 

{^h)22 ~ ^ tp'i {^h)is — ^v{2fiefiSi3 — (Aa -|- (13) 

{^11)23 ~ ^v{2p,esCj3 — {Ax + p!)sp)] = (Ai;^(Aa -|- p!) — {^s + 
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where fieS = As, tjs = tan/3 and r = /iefr(^A + 3^0 + The CP-odd Higgs mass 

matrix in the basis is given by 


(^l)ii = 2r/s2/3; {Ma)i2 = 

(^^ 1)22 = s (A'1'^(-4 a -h m')'S/3C/9 - + ^s)) - ‘^rn'l, 


(14) 


and the charged Higgs mass is 


Ml = 2r/s2y3 - (A^ - gl/2)v\ 


(15) 


We point out two features of the tree level masses that will be important in our 
discussion of the one-loop corrected Higgs mass. The first feature is the correlation 


among the scalar masses in the limit where Ax and are small compared to the heavy 


Higgs masses. This is best seen by setting (M^) = 0 in Eq. (14) (which can be obtained 



(dH), we find that the CP-even, CP-odd and charged Higgs eigenstates arising from the 
SU{2) doublet sector are nearly degenerate in mass, a feature well-known in the MSSM. 
Their mass splittings ~ These three fields then have a mass ~ Ma^^ in the limit 
/§> where denotes the corresponding CP-odd eigenmass. Likewise, 

the CP-even and CP-odd Higgs eigenstates arising from the SU{2) singlet are nearly 
degenerate, with mass splitting ~ s^. Therefore, these two fields have a mass ~ MAg in 
the limit M^^ /§> 

The second feature is the decoupling of heavy states. Raising M^^ and M^^ decreases 
their impact on the mass of the lightest CP-even state, effectively making it more SM- 
like. A simple way to see this decoupling behavior is to rotate the CP-even mass matrix 
into the basis 


= hls/s -h = hlcp - h^s = K 


(16) 


which leads to the CP-even mass matrix 


{^H)hh ~ i^H)hH ~ (A^'i’^ - -Tf|)s4^/2; 

{^^h)hh ~ ^Ad ~ (A^3^^ — AL|)s|^; = 2Xv{fieS — AaS2/3); 


(17) 



Notice that and are absorbed into our definition of M^^ and M^^. For large 


tan/3, MI^ and is identified with the SM Higgs, and and with non¬ 


standard Higgs bosons. This decoupling feature should be preserved after the inclusion 


of radiative corrections to the lightest CP-even Higgs mass, which is a non-trivial check 
of this computation. 
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2.2.2 Radiative corrections 


The mass of the lightest Higgs boson can be significantly modified by one-loop correc¬ 
tions. The largest contributions to the Higgs potential at one-loop level are from the 
Higgs bosons, third generation squarks, charginos and neutralinos. Thus we have 


AH 


1 

327r2 




6AH‘ 





where for the sector in the MS scheme, 

i ^ ' 2 =heavy 



.(19) 


Mfa{4>k) is the field-dependent mass eigenvalue for contribution, Mf^ is the corre¬ 
sponding field-independent tree-level eigenvalue and the renormalization scale Q ~ mh = 


125 GeV. The approximation in Eq. (19) holds because we are interested in large correc¬ 


tions to the lightest Higgs mass due to states much heavier than the electroweak scale. 
Also, the field dependences inside log terms are neglected since they only induce higher 
order field-dependent terms. 

The dominant contributions to AH in our scenario are due to heavy Higgs scalars 
coupling to the light Higgs boson with strengths proportional to powers of A. The effects 
of the top quarks and the scalar tops on the Higgs potential have been studied in great 
detail in Refs. [H El El d El [U EDI EH EH ESI ESj- To highlight the effect of large A 
corrections, we suppress the contribution of scalar tops to AH by choosing their masses 
close to the electroweak scale while still being compatible with ATLAS and CMS bounds 
HZIISSIISS]. The contributions of charginos and neutralinos to AH are typically small. 
The Higgs couples to the bino and the wino triplet with electroweak strength, whereas 
the A-dependent coupling to the Higgsinos and singlino is typically suppressed due to 
neutralino mixing. In addition, the masses of the Higgsinos and singlino ^ 1 TeV in 
our phenomenology while Ma G [4,8] TeV. We therefore neglect corrections from the 
chargino-neutralino sector in the remainder of this article. 


In order to compute the one-loop corrections to the Higgs potential in Eq. (19) due 


to the heavy CP-even, CP-odd and charged Higgs bosons, we must determine the field- 
dependent eigenvalues of each of the respective matrices. When expressed in terms of the 
matrix elements these field-dependent eigenvalues can in general be quite complicated. 
The calculation can nevertheless be simplified if we expand the eigenvalues as a Taylor 
series in the tree-level masses: 


Mi{4>k) — Mftj-ee + bi{4>k) + 

-IVlj , 



( 20 ) 


where the coefficients 6i(</>fc) and Ci{(pk) are at most quadratic and quartic in the fields 
respectively. Purthermore, when evaluated at the tree-level VEVs, the coefficients satisfy 
the condition b{{vk}) = 0 = c{{vk})- 

In practice, we expand the eigenvalues as a Taylor series in the pseudoscalar masses 
and Using these approximations the one-loop effective potential due to the 









heavy Higgs scalars now evaluates to 


AH oc ^ [a^MXi + 2biMl, + (bf + 2q)] log 


m2 


A,i 




( 21 ) 


where a* are constants and the field-dependent coefficients bi and c* are obtained from 


the hatted coefficients in Eq. (20). Reducing AH to this form considerably simplifies the 


calculation of Higgs mass corrections. AH as presented here must also ensure that the 
decoupling behavior discussed in the previous subsection is preserved at one-loop order. 
This result is demonstrated in Appendix 

The full details of our computation and the corresponding results are presented in 


Appendix where two cases satisfying the condition (M^) = 0 in Eq. (14) were 

considered. In the hrst case, which we call Case (A), we assume that the one-loop 
corrections arise from a single heavy scale Ma = M^^ = M^g- The results from 
this case will be used in our discussion of phenomenology in Section In the second 
case, which we call Case (B), we show the effect of splitting the CP-odd Higgs masses, 
thereby obtaining corrections from two heavy scales. In this case we set the terms 
Ax, fj,', m 3 , m'g to zero for simplicity. Further, we ignored corrections that depend on 

electroweak couplings since we are interested in the limit g. It is important to note 
that Cases (A) and (B) pertain not only to different limits of the mass spectra of the 
CP-odd scalars, but also to somewhat different regions of the Lagrangian parameters. 
In Case (A), the parameters Ax, g!, Ag^,m^,m'g can be non-zero in general, with the 
condition = 0 imposing Ax = g'■ On the other hand, Case (B) explicitly sets 

them all to zero. 

For Case B, the one-loop self-energy corrections obtained in the basis (h)), hj), /i°) are 


Hii = 


ni2 = 


(Ml \ [Ml' 

-( 4 c 2/3 -h C 4/3 + 1 ) log ( 1 -h 2 log ' ® 


AS2 2 
167r2 ^ ^ 

2AS2 (Ml 


Af| 


Ml 


„ 2 

1122 — - 

167r2 ^ 

A^U/Xeff 

IGtt^ 

A^U/Xeff 

IGvr^ 

4A% 

167r2 


-(-4c2/3 + + l)log 


'yh' 

Ml 


-k 21og 


His = 


n23 = 


Has = 


sp 


Cfi 


-(1 - 3C2/3) log 


Ad 

M| y 

1 -k41og 

Ml \ 


Ml 1 

1 -k41og 



Af| 


log 


'yh' 

Ml 


When these contributions are rotated into the basis of Eq. (IG), we get the (1, 1) element 


9 






























\ / 
h \ / h 

\ / 

\ / 

X 

/ \ 

/ \ 


^heavy^ 

/ 


\ 

\ 


(a) 


/h 

X K 

h / ,*!“'■>’ ^ h 

/ * \ 

(b) 


Figure 2: (a): Tree level quartic vertices involving at least two h fields with vertex factors oc A^, 
in the limit tan/3 3> 1. In this limit, —>■ /i, —>■ H, h^. No quartic terms at formed at 

tree level. (j)i correspond to the heavy fields H,h^, Ag. (b): One-loop quartic vertices with 
four h legs, formed from the tree level vertices in (a). These are oc and account for most of 
the radiative corrections to the Higgs mass in our model. 


of the self-energy corrections as 
X^v^Sjs 


^hh = 


IGtt^ 


f‘% 


(c^(2 + C 4 ^) — S^(l -I- C 4 ^ -h 4 c 2/3)) log ( j -I- 2s^ log 


Ml 


( 22 ) 


This is a good approximation for the Higgs mass correction when the mixing between 


the SU{2) Higgs doublets and the singlet is negligible. At large tan/3, Eq. (22) further 
simplifies to 


large tan/3 


(Ml \ [Ml' 

21ogK# +21og' 


Ml 


Ml 


(23) 


We could gain an intuitive understanding of our results by qualitatively estimating 
the size of the one-loop radiative corrections without recourse to the effective potential. 
Such an estimate would serve as a useful cross-check of the results obtained from AF. 
We do this by the following argument in our limit of interest, tan/3 /!> 1 and Xl$> g. In 
this limit, we identify the real scalars —?• /i, ^ H, —?• h^, where h is the SM 
Higgs boson, and H and are non-standard Higgses. The Standard Model Higgs and 
the Goldstone bosons reside mostly in and the non-standard CP-even and CP-odd 
Higgses in and S. 

For X g, the most important quartic terms at tree-level are those proportional 
to A^. Before EWSB, we can read them off from Eq. (10) as the terms 
F0|2|5|2 and 




\HT 


I tr 0 | 2 | 17012 


|5p. After EWSB, we can expand Hu,Hd,S using 
Eq. © to obtain various quartic vertices in terms of the real and charged scalars. 

Fig. j^a) shows all the tree-level quartic vertices that involve at least two h fields. 
Recall that the SM Higgs mass is set by the coupling strength of the quartic term h‘^ in 
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the scalar potential. The tree-level A-dependent quartic terms are suppressed at large 
tan/?. However, using the vertices in Fig. [^a), we can construct four one-loop level 
quartic vertices proportional to as shown in Fig. lib). Each of these diagrams is 
proportional to log(M^,/M^), where is the mass scale of the heavy field running 
in the loop. Two diagrams each correspond to and respectively. Since the 

internal propagators are identical, each diagram comes with a factor of 2. Canonical 
normalization of the mass term of a real scalar implies an additional factor of 1/2. 
Finally, including the loop factor we find the approximate correction to the 

lightest CP-even eigenstate to be 


n 


hh 


1 2 • 2 • A^ 

2 ’ 167r2 


log 


'Ml \ fMl ' 
^ol+logl "Is 


M| 


Ml 


(24) 


which agrees with Eq. (23). 


It would be interesting to compare the Higgs mass corrections obtained from the 
heavy Higgs fields and those obtained from top squarks. For simplicity, let us set the 
pseudoscalar masses equal, Ma = = Maq , and obtain 


T-rhiggs _ 


47r2 


log 


yii 

Ml)' 


(25) 


Again for simplicity, we can assume the top squarks are degenerate (m^ = = mf^). 

Then we obtain [33] 


H 


stops 

hh 




(26) 


le pseudoscalars and the top 

andHlth«tnJf'>nfr 


The factor of 3 arises from the three QCD colors. If t 
squark are degenerate {Ma = from Eqs. (25 

for A > 3^/^yt. Since yt = mt/v ~ 1, we have for A > 1.3 

In the discussion of our model’s phenomenology, we set = 800 GeV while Ma 
ranges between 4 TeV and 8 TeV; therefore, the one-loop corrections from the Higgs 
sector dominate those from the stops. Hence, throughout our analysis, the effect of the 
top squark correction to the SM Higgs mass is neglected. 

We can now quantify the discrepancies between the results obtained by a full one- 
loop effective potential calculation and those obtained by other means. To do so, first we 


compute the correction to the Higgs squared mass obtained from Eq. (22), and denote 
it by Am^. Eor the same set of parameters, we compute (Am^)j for each alternative 
approximation labelled by i. We then take the difference and normalize it to Am| and 
define the discrepancy as 


( 5 (Am|) = 


(Am^)j - Am^ 


Am| 


(27) 


which is then expressed as a percentage. This approach eliminates the A-dependence of 
the discrepancies and allows us to focus on their behavior with respect to tan /? and the 
heavy (pseudo)scalar masses. 
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Figure 3: Discrepancies between the Higss mass radiative corrections obtained from our one- 
loop effective potential in Eq. (221 and those obtained by other means, as a function of the mass 
Ma of degenerate pseudoscalars. The blue, dashed red and magenta curves represent corrections 
obtained from Eq. (221, Eq. (24) and Ref. [3S] respectively, (a) corresponds to tan/3 = 2, (b) 
corresponds to ia,nj3 = 50. See text for details of the behavior of the curves. 


Assuming for simplicity that the CP-odd scalars are degenerate, we depict in Fig. 
the discrepancies as a function of Ma- Figs, [^a) and[^b) correspond to tan/3 = 2 and 
tan/3 = 50 respectively. The blue curve denotes (Am|)j obtained from the approxima¬ 
tion in Eq. (22). Since this approximation neglects doublet-singlet mixing, it tends to 


overestimate the correction, i.e., (5(Am|) > 0 as observed in the plot. The discrepancy is 
also seen to asymptote to zero at large Ma-, where the CP-even singlet Higgs decouples 
from the SM Higgs. The dashed red curve is (Am^)j obtained from our qualitative dia¬ 


grammatic estimate (Eq. (23)). Since the estimate is designed for large tan/3 it disagrees 
with the blue curve at tan/3 = 2, but coincides with it very well at tan/3 = 50. The 
magenta curve depicts (Am^)j obtained from NMSSMTools 4.5.1 [35l 03], which also 
computes the one-loop radiative corrections from the effective potential, albeit under a 
different set of approximations. We find an interesting discrepancy here, to which we 
now turn. 


The eigenvalues of the CP-odd mass matrix in Eq. (14) are given by 


E‘i = - (T±y/W^^ 


(28) 


where T = -|- (M ^)22 trace and D = (M ^)22 “ {^a)‘i 2 

determinant of the mass matrix. In Ref. [35], it is assumed that D T^, so that 
the eigenmasses are obtained as E'^ ~ T, E'^ D/T. This always leads to a hierarchy 
between the pseudoscalar masses. In contrast, our approach sets the off-diagonal element 
(M ^)^2 to so that the eigenmasses are = M^^ = 

(M^) 22 - Therefore, our approach allows for a variety of mass splittings. Hence the 
discrepancy between us and Ref. m is expected to be maximum when the CP-odd 
Higgses are degenerate, and minimum when these masses are well split. We illustrate 
this effect in Fig. Since (M ^)^2 = 0 in our approach, we set to ^o^^o in 

the expression of Ref. |3S] as well, in order to make an “apples-to-apples” comparison. 
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Figure 4: The discrepancy between Higgs mass corrections obtained by Ref. |35] (which were 
used in the original code of NMSSMTools 4.5.1) and by us, as a function of the ratio of the heavy 
CP-odd Higgs masses. The red (blue) curve corresponds to M^s = 1(2) TeV. The discrepancy 
arises due to an approximation assumed by Ref. |35j . namely, that a hierarchy exists in the 
pseudoscalar spectrum. It is seen that our results agree when there is indeed a hierarchy. See 
text for more details. 


= 1 TeV 
Ma, = 2 TeV 


We then plot as a function of Ma^/^As^ where we have taken A = 1.25, 

tan/3 = 50 and fieS = HO GeV. The red and blue curves depict Mas = 1 TeV and 
= 2 TeV respectively. As expected, we find the discrepancy at its greatest at 
Ma^/Mas = 1, which can reach upto ~ 15%. Observe also that (5(Am|) < 0, implying 
that Ref. |35j underestimates the one-loop contribution to the Higgs mass in the region 
around Maj^/Maq = 1- As we raise the discrepancy drops quickly and our 

results concur. 

The results of Ref. [33] were originally used in the code of NMSSMTools 4.5.1 [33] , 
Since our phenomenology in Section assumes Ma^ = Mas , we replaced the code in 
NMSSMTools 4.5.1 with the expressions that we derived in Appendix |B| 

2.2.3 Stability of the electroweak scale 

The minimization conditions of the tree level Higgs potential in Eq. ([^ lead to the same 
relation between the electroweak scale and the SUSY parameters seen in the MSSM. In 
particular, the EWSB condition is |44j 

^2 _i_ 1 

= ,2 _ -I ^ iFeffl^ ) (29) 

1/3 i 

which at large tan j3 reduces to 

-M| Ri — Ih-effl^ , (30) 

where the terms are suppressed by With this result we can now quantify the 
relative importance of different contributions (denoted by a) to the EWSB scale (M^/2) 
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as 


A(a2) 



(31) 


The tree-level and one-loop corrections are the same as in the MSSM and are well- 
known [l5]. For instance, the tree-level contribution due to 350 GeV is equivalent 

to the one-loop contribution of stops at ^ 800 GeV |l 6 ] . Hence the regions we are 
considering in this article are typically as tuned as regions of the MSSM with a light 
stop. 


2.2.4 Higgs couplings to SM particles 

LHC measurements of signal strengths (production rate x branching ratio) can poten¬ 
tially constrain the properties of the Higgs sector. Mixing among the Higgs fields can 
in principle alter the lightest Higgs boson’s SM-like behavior. We follow the analysis of 
Ref. [30] to apply the relevant limits. 

After including the one-loop self-energy corrections, we rotate the Higgs fields (/i°, h^, h^) 
into the mass eigenbasis (/ii, /i 2 ) ^ 3 ) and identify the lightest scalar as 

hi = {—h^ sina + h^ cos a) cos 7 -|-/i^ sin 7 , (32) 

where the angles a is the usual MSSM CP-even mixing angle that characterizes doublet- 
doublet mixing and 7 characterizes the doublet-singlet mixing. We can then write down 
the reduced couplings of hi to pairs of fermions and vector bosons as 

9tth\ 

= cos 7 (cos 5 — sin (5 tan/3), 

dbbhsM 

9vvhi f- / 0 Q\ 

- = cos 7 COSO, (33) 

9VVhsM 

where 6 = a — P + 7r/2. 

If we inspect the off-diagonal entries of Eq. dlT] ), we see that for Ax <C Ma and large 
tan/3, (M^)^^ < (M|^)^^. Thus as we raise Ma, the heavy doublet Higgs (identified 
as /13) generally decouples faster than the heavy singlet (identified as / 12 ), as noted by 
Refs. [251 El] In dealing with the phenomenological consequences of our model, we focus 
exactly on the region of Ax <C Ma and large tan/3. Therefore for the rest of this analysis 
we assume h^ is decoupled from the spectrum and /12 is not. In this limit, the mixing 
angle 7 is given by 


cos 7 ( cos S + 


sin (3 
tan/3 


. 2 

sm 7 = — K — „ 
^2 - 


hi 


(34) 


where sin^ 2/3 -|- cos^ 2/3, and the Higgs couplings to fermions and vector 

bosons become 


gtthi 9bbhi gvvhi 

-=-=-= cos 7 . 

dtthgM dbbhQyi 


(35) 
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Using these relations Ref. m performed a universal fit on the LHC signal strength 
measurements and found that sin^ 7 < 0.23 at 95% C.L. This result was obtained using 
tree level relations for the reduced couplings. When we include our one-loop corrections, 
we find that the reduced couplings are modified by less than 1%. Therefore, in the 
discussion of our model’s phenomenology in Section we will simply use the results of 
Ref. |3U| to constrain the Higgs couplings with LHC measurements. 


2.3 Neutralino sector 


The composition of the lightest neutralino and its couplings to the Higgs sector is central 
to the dark matter phenomenology of our model. The neutralino mass matrix in the basis 
{B, W, is given by 


Atneut — 


/ Ml 0 —giv cos P/Vl givsin/3/V2 0 

0 M 2 g2VCOsj3/y/2 — g2V s\n j3/y/2 0 

—cos/3/\/2 g2Vcos/3/y/2 0 —geS —Ausin/3 

5 iusin/ 3/\/2 —g2Vsinl3/y/2 —g^eS 0 —Xvcos/3 

\ 0 0 — Ausin/3 —Xvcos/3 g' 


Notice that when g' <C Mi, M 2 , ges, large A couplings imply a large Higgsino component 
in the lightest neutralino. This feature has many unique consequences for the dark matter 
phenomenology discussed in Sec. As we shall see, the Higgs-Xi-y? coupling strengh 
plays an important role in constraining our model with dark matter experiments. This 
coupling, denoted hereafter by ghxxi obtained as 


A 


9hxx = 






(37) 


where the W and (/j are the appropriate components of the lightest neutralino and the 
SM-like Higgs respectively. In terms of the rotation angles in Eq. (32), we can read off 


= — sinacosy, = cosacosy, Cs = sin 7 . 

The dominant channel for Xj’-nucleon scattering is through a t-channel Higgs. There¬ 
fore, dark matter direct detection experiments, as well as limits on the invisible decay 
width of the Higgs, apply strong contraints on the coupling ghxx- ^ suppressed ghxx 
occur in our model either when the Higgsino content is suppressed, making Xi niostly 
singlino or bino, or when there is a delicate cancellation between the various terms in 
Eq. 1^ We illustrate this point in more detail in Sec. 


2.3.1 Electroweak precision limits 

Due to mixing between the Higgsinos and the singlino induced by large A in certain 
regions, constraints from electroweak precision experiments can be strong in Fat Higgs/A- 
SUSY models [221124| . In particular, the T parameter can get large contributions from 
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Figure 5: Limits from electroweak precision parameter T on the neutralino sector of our model. 
The shaded regions are where > 0.15 and therefore excluded at 95% C.L. Regions shaded gray 
correspond to the wino decoupled from the spectrum (M^ = 10 TeV) and regions shaded red to 
= 200 GeV. In (a), A = 1.25 and tan/3 = 5 and in (b), /Xes = fi' = 300 GeV. See text for 
details of the behavior of these curves. 


the neutralino sector, denoted hereafter by T^. This phenomenon is understood easily 
in the limit where the electroweak gauginos B and W decouple from the spectrum, i.e.. 
Ml, M 2 are very large. This leaves us with three mass scales and \v, which set 

the mass of the lightest neutralino, M^o. The lightest chargino is mostly Higgsino with 

Al 


a mass /ieff- In this limit, is large when Tf-± — M^o is large and when there is as 




a significant Higgsino component in Xi- For simplicity, let us work in the limit where 
tan/3 is large. Then the neutralino mass matrix in Eq. (36) is simply 


A^neut ~ 


^ 0 -/Teff -Ar; 

/^efF 0 0 

—Xv 0 


(38) 


Ty, is suppressed either when /r' ~ /^eff ~ Au, where M^± — M^o is small, or when 
fJ-eS ^ ^ Xv, where the Higgsino component in Xi is suppressed. For //' <C HeS ~ Au, 

where both M-± — M^o and the Higgsino component in v? are large, constraints from 
Ty. can be strong. 

Lowering the mass of the wino triplet M 2 to ~ /^eff ~ Au can have a significant impact 
on Tjj.. This is because the wino would mix with the light neutralinos and charginos. 
Lowering the bino mass Mi, on the other hand, gives only a negligible contribution to 
Ty-. This is because the bino mixing with the rest of the neutralinos is only proportional 
to 510 

^ It must be remembered that relative minus signs between and would introduce quantitative 

changes in the picture owing to new phases in the neutralino mixing angles. We will not include these relative 
signs in our discussion. 
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In Fig. we present the T-parameter contributions from the charginos and neu- 
tralinos, which were computed using the general expressions provided in Ref. |50j . In 
Fig. I^a), we take A = 1.25 and tan/I = 5 and show our results in the /Teff — /x' plane. 
The shaded regions denote where T^- is not within the 95% C.L range [—.01,0.15] set 
by the Particle Data Group m- The gray region corresponds to large gaugino masses 
{Ml, M 2 ) = (10 TeV, 10 TeV) while the red region corresponds to a light wino with 
{Ml, M 2 ) = (10 TeV, 0.2 TeV). As discussed above, lowering the wino mass can lead to 
a larger T^. For small //', decreases as /Xeff increases due to a reduction in the Higgsino 
component of the lightest neutralino. Similarly, raising /x' has the effect of reducing the 
splittings between the neutralinos and charginos which also leads to a smaller T^. 

The effects of varying A and tan/3 on are presented in Fig. [^b). Here we fix 
/Xeff = /x' = 300 GeV. The colored regions have the same definition as those in Fig. [^a). 
Since the elements of A^neut quickly asymptote to fixed values as a function of tan /3, it 
can be seen that is insensitive to large tan /3. This insensitivity to large tan /3 is clear 
in the relation derived in Ref. |22j 


Tx 


(izl 

U+1 


2 

A(/Xeff , fJ- , A), 


(39) 


where T(/Xefr,/x', A) is some function of these variables. This relation also shows that 

is suppressed as tjs approaches 1, thereby allowing for larger values of A. As stated 
before, lowering M 2 typically increases the neutralino and chargino contributions to the 
T-parameter. However, it is important to emphasize that increasing either /Xefr or /x' can 
significantly lower the electroweak precision constraints even for large tan/3. A large ;Ueff 
comes at the cost of a slight increase in electroweak fine-tuning, but can greatly weaken 
T-parameter constraints. 

Finally, we make two remarks. First, the 5-parameter was not discussed here. This 
is because the contributions of our model to 5 are very small in our regions of interest 
and hence the constraints are much weaker than those on the T-parameter. Second, the 
T-parameter receives a stop-sbottom contribution, as discussed in Ref. [22]. In the limit 
of zero left-right mixing, this is given by 


Tst-sh 


500 GeV\ 

J 


(40) 


In our phenomenological discussions, we will choose = 800 GeV to suppress this 
contribution. 


3 Phenomenology 

In this section we study the phenomenological constraints on the large tan/3 regions 
of the Fat Higgs/A-SUSY models. In addition to the constraints arising from Higgs 
corrections discussed in the previous section, we also include limits from dark matter 
experiments, most importantly those set by the LUX experiment |52] . In particular, 
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the mass and couplings of the lightest neutralino Xi can put strong constraints on our 
parameter space. 

In order find phenomenologically viable regions, we modihed NMSSMTools 4.5.1 [l3] 
to include the Higgs mass corrections we computed in Sec. 2.2.2 We then made the 
following simplifying assumptions; 

• In the Higgs sector, we take the pseudoscalars to be degenerate, with = 

M^g = Ma- Furthermore we assume that m'^ = m 3 = 0, so that the heavy CP- 
even Higgs bosons are also (nearly) degenerate. The condition that the CP-odd 
masses are degenerate requires that (M^) = 0 in Eq. 


14 


which implies Ax = /ih 


Therefore, both /r' and fieS control the amount of doublet-singlet mixing in Eq. (13). 
The only independent parameters in the Higgs sector are then: A, fieS, /u',tan/3 and 
Ma. 

In order to be safe from electroweak precision bounds, we decouple the winos at 
M 2 = 10 TeV, leading to an effective theory for the neutralino system with five free 
parameters : Mi, //eff, A and tan/3. 


We require /ieff > 104 GeV to evade the LEP H bound on charged Higgsinos 

The sleptons and the first two generations of squarks are decoupled from the low 
energy phenomenology and their masses set at 5 TeV, unless stated otherwise. The 
top squark parameters are set at m^^ = = 800 GeV and At = 0, thereby 

making the stop contributions to the Higgs mass and the electroweak symmetry 


breaking condition in Eq. (30) small. This choice of stop masses also avoids con¬ 
straints from collider searches @711181119] and, as mentioned in Sec. |2.3.11 from 
electroweak precision tests. 

• We choose to require the conventional upper limit tan/3 < 60, so that 2/6 < 1 at 
the weak scale. Larger values of yt, may be allowed as long they do not develop a 
Landau pole at a scale below that of A. 

• We assume that the relic density of dark matter is the value determined by Planck 
|53j . Hence, in scenarios where < 0.12, we assume some nonthermal mecha¬ 
nism for generating the observed relic abundance. 

These assumptions reduce the number of independent SUSY parameters to 

A, tan/3, Ma, fJ-eS, IJ-', Mi . 


As discussed in Sec. 2.3.1, precision electroweak constraints are weak either when /r' or 
/ieff are large for any tan/3, or when /i' <C Au ~ /ieS at low tan/3. In these regions, ghxx 


(as defined in Eq. 37) can also be found to satisfy dark matter direct detection and relic 


density constraints. In particular we find the following viable parametric regions. 


3.1 Singlino DM: Large tan^^ and /i' < At ~ /igff <C Mi 

Large A and large tan j3 are a new region of parameter space that have not been empha¬ 


sized in the literature before. We showed in Sec. 2.2 that this region can be compatible 
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with the mass of the SM Higgs boson because one-loop radiative corrections to the Higgs 
mass are insensitive to tan/? at large values, and are set solely by A and Ma- We also 
showed that precision electroweak constraints can be weak in this region. We now show 
that this region is also compatible with constraints from dark matter. 

As mentioned in Sec. 


2.3 


it can be seen from Eq. (37) that Qhxx suppressed when Xi 
is mostly singlino such that N~^ ~ 1. This requirement is possible when fi' is relatively 

‘^s 

small compared to the other mass scales in the neutralino mass matrix. The annihilation 
of Xi into SM fields in the early universe is generally inefficient, due to both the Z- and h- 
mediated channels being suppressed by the small Higgsino component of Xi- Therefore, 
for the cosmological relic abundance to be below the observed value ~ 0.12, we 

consider the mechanism of resonant annihilation and co-annihilation 1591. 


3.1.1 Resonant annihilation 

If M^o happens to be close to m/i/2, it can undergo resonant annihilation through an 
s-channel Higgs. Therefore, we set = 62.5 GeV in this scenario. We also set /Xeff = 
800 GeV, Ml = 1 TeV, tan /? = 50 and Ma = 4 TeV. The orange curves in Fig. j^a) 
depict contours of the LSP-nucleon scattering rates, usi (in units of uo = 10“^® cm^), 
on the fieS — A plane. The red shaded regions are excluded by LUX at 90% G.L., and 
the green band corresponds to 120 GeV < rrih < 130 GeV. Contours of are denoted 
by dashed curves. 

In Fig. [^a), the dark matter-nucleon scattering rates are smaller for larger /igfr be¬ 
cause the Higgsino fraction in Xi decreases. is observed to rise with increasing A due 
to an increase in the Higgsino fraction of The region around rrih ~ 125 GeV corre¬ 

sponds to ~ 0.05, which is safe from electroweak precision constraints. This regions 
is also safe from invisible Higgs decay bounds since the process h ^ XiXi is phase space 
suppressed. 


3.1.2 Co-annihilation region 

If the mass spectrum is such that one or more sfermions are nearly degenerate with 
Xi, dark matter annihilation could be assisted by the sfermions through co-annihilation 
effects, leading to a small relic abundance. Bounds from LEP on charged sfermions m 
would then imply that M^o ^104 GeV. 

We investigate this scenario in Fig. [^b). In this figure we assume that the correct 
thermal dark matter relic abundance is generated by a process like co-annihilation. We 
do not explicitly state the mass spectrum or compute the resultant relic abundance. 
Again the orange curves in Fig. |^b) are contours of crgi (in units of do = 10“^^ cm^) 
on the fieS — plane. We have chosen A = 1.25, tan/3 = 50, Mi = 1 TeV and Ma = 8 
TeV. The larger values of /Xefr compared to those in Fig. [^a) imply a greater amount 
doublet-singlet mixing. Therefore a slightly larger value of Ma is chosen in this scenario 
as compared to that of Fig. |^a). The region shaded red is excluded by LUX at 90% 
C.L. The green (gray) vertical bands correspond to rrih G [120,130] GeV for A = 1.1 
(1.25). The effect of varying A on dsi is not shown since the scattering cross-section is 
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Figure 6: Large tan/3 parametric scenarios for Fat Higgs/A SUSY models, (a) Singlino resonant 
annihilation {p! = 62.5 GeV, Mi = 1 TeV): orange curves correspond to crsi in units of (Tq = 
10“^^ cm^. The red shaded region is excluded by LUX at 90% C.L. and the green shaded region 
corresponds to 120 GeV < mh < 130 GeV. The dashed lines are contours of T^. (b) Singlino 
co-annihilation: orange curves and red region the same as (a). The green and gray shaded regions 
correspond to 120 GeV < mh < 130 GeV for A = 1.1 and A = 1.25 respectively, (c) Bino resonant 
annihilation (/r' = 1 TeV, Mi = 62.5 GeV): blue curves correspond to ctsi, with the remaining 
colors the same as in (a), (d) Bino co-annihilation: blue curves correspond to crsi, with the 
remaining colors remaining the same as in (b). In all these plots we have set tan/3 = 50. The 
critical features are explained in the text. 


20 



















insensitive to it due to the large values of /Xeg. The decrease in usi with /Xeff is due to 
the decoupling of the Higgsinos leading to the suppression of ghxx- increase in usi 
with g! is due to the larger Higgsino fraction in Xi, which leads to an enhanced Qhxx- 
The relatively large size of and g! here suppress the Higgsino sector contributions to 
the T parameter. For regions where m/j ~ 125 GeV, we find that < 0.03. 


3.2 Bino dark matter 


This parametric scenario is the bino analogue of the previous singlino dark matter sce¬ 
nario we have discussed. It is also a typical scenario that arises in the MSSM at large 
tan/?, ghxx-i ™ ^0- (37), will again be suppressed when Xi is mostly bino. The bino 
fraction of Xi is increased by lowering Mi relative to other mass scales in the neutralino 
mass matrix. Again the observed cosmological relic abundance is either through the 
mechanisms of resonant annihilation and co-annihilation 1591. 


3.2.1 Resonant annihilation region 

Again in the resonant annihilation region, we set Mi = 62.5 GeV, //efr = 800 GeV, 
fi' = I TeV, tan/3 = 50 and Ma = 4 TeV. The blue curves in Fig. [^c ) depict contours 
of the LSP-nucleon scattering rates, usi (in units of do = 10“^® cm^). The remaining 
colored contours correspond to the same regions as those in Fig. |^a). The scattering 
cross-sections are stronger here than in Fig. [^a) because in the singlino-like scenario 


there is a partial cancellation among the terms in Eq. (37), which suppresses ghxx- This 
cancellation arises from an extra minus sign picked up by N^q for the range of mass 
parameters considered. 

Similar to Fig. [^a), the dark matter-nucleon scattering rates are seen to decrease 
as we decouple the Higgsinos by increasing /Xefr- In contrast to singlino dark matter, 
Tjj, ~ 0 for bino dark matter throughout the plot in Fig. j^c) because both the charged 
and neutral Higgsinos are quite degenerate. 


3.2.2 Co-annihilation region 

Similar to the singlino scenario, we assume that the sfermion mass spectrum is such 
the relic density of x? is consistent with cosmological observations. The blue curves in 
Fig. [^d) depict contours of dsi in units of do = 10“'^® cm^. We vary Mi while fixing 
/r' = 1 TeV, and the remaining parameters are the same as in Fig. j^b). The dependence 
of dsi on |Ueff and Mi is similar to that of singlino scenario with g,' —)> Mi. Since the 
Higgsino fraction is larger in the mostly singlino x? that the mostly bino Xi^ is large 
in Fig. I^b) compared to Fig. [^d). For regions where m/j ~ 125GeV, we find Ty, ~ 0. 
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A=0.9, tanj8=2.5 


A=0.75, tan/3=1.5 
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Figure 7: The well-tempered scenario at low tan/3, with Xi an admixture of bino, Higgsino 
and singlino. In (a), A = 0.75, tan/3 = 1.5 and in (b), A = 0.9, tan,d = 2.5. The heavy Higgs 
states are decoupled at Ma = 5 TeV. This choice of parameters fixes ~ 125 GeV. Regions 
shaded red are excluded by LUX at 90% C.L., blue hy h ^ Xi Xi bounds and gray by Z —>• y? Xi 
bounds. These constraints leave a small patch of parameter space that are still viable, the “blind 
spots”. The dashed lines are contours of M^o in GeV. More details are presented in the text. 


3.3 The well-tempered bino/singlino/Higgsino 

In the limit where fj,' <C fieS, Mi, precision electroweak contraints can be evaded by 
raising thereby decoupling the Higgsinos. However, raising or tan (3 suppresses 
the mass x\ a-s 

M^o Ri /i' -b A^u2/refrS2/3/(/Ueff + (41) 

for large Mi and M 2 . As M^o < Mzj^ for a large region of parameters in this scenario, 
the invisible width of the Z boson is an important constraint. Consequently, to find a 
viable region of parameter space, we require /iefr ~ An and small tan/3. In this region, 

Ohxx supppressed when y? is an admixture of B, and ipg such that they lead 
to “blind spots” in parameter space [12] - regions that are compatible with current 
experiment. For illustration, we have consistently set /u' = 0 in this section. 

We illustrate these blind spots in Fig. which shows constraints on the LSP in the 
Ml — /ieff plane. Fig. j^a) corresponds to tan/3 = 1.5 and Fig. Ob) to tan/3 = 2.5. 

To fix mh ~ 125 GeV, we take A = 0.75 and A = 0.9 respectively and decouple the 
heavy Higgses with Ma = 5 TeV. At these values of A, tan/3 and Ma the Higgs mass is 
mainly set by the tree-level values as the loop level corrections are small. The regions 
shaded red are excluded by LUX at 90% C.L. Regions shaded blue are excluded by the 
latest limit on the invisible decay of the Higgs, B.Tl.{h —>• XiXi) < 0-44 (57] [58|. The 
gray region is excluded by limits from the invisible decay of the Z. The dashed curves 

^As an aside, any admixture of bino and singlino such that A~ -|- ~ 1 will also lead to a suppressed 

9hxx- 
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represent contours of M^o in GeV. This range of parameters is cosmologically viable with 

A1 

< 0.12, where the dominant primordial annihilation of Xi is through an s-channel 

Z. 

A comparison across the plots informs us that an increase in tan /3 strengthens the 
constraints from Z, h —>■ XiXi) which is due to the decrease in M^o, as discussed earlier. 
We also notice that the LUX constraints are consistently stronger than h —)• XiXi bounds. 
Therefore, the blind spots (unshaded regions) are determined in this case by limits 
from LUX and invisible Z decays alone. As mentioned in Sec. |2.3[ larger values of A 
contribute more to T^. For Fig. [^a) and Fig. [^b), < 0.02 (completely safe) and 

< 0.11 (marginally safe) in the blind spots. Also, using NMSSMTools 4.5.1 we find 
that ~ 0.01 in the blind spots. Therefore, the well-tempered neutralino here can 

only make up a small fraction of the observed relic density. 

3.4 Future prospects 

In this region (large tan f3 with ~ Au ~ figs ), the non-standard Higgs scalars are heavy 
with Ma between 4 — 8 TeV. Therefore the doublet-singlet mixing in the Higgs sector is 
very small leading to a 1% deviation in the Higgs signal strengths from the SM. Such 
deviations are much below the sensitivity of the LHC at present and future runs, and 
can only be tested at a future “Higgs factory”. However, the large tan (5 scenarios can be 
probed by future dark matter direct detection experiments. In particular, the projected 
reach of the XENONIT experiment |36j corresponds to usi ~ 10“'^^ — 10“^® cm^ for 
dark matter masses between 50 GeV and 500 GeV. Since the DM-nucleon scattering 
cross-sections for the large tan /3 scenarios in Fig. vary from rsj 10-46 _ 10-45 cm2, 
these regions can be probed at the XENONIT experiment. Unlike the large ^',tan/3 
scenarios, XENONIT will only be able to probe some of the allowed regions of the 
well-tempered scenario because ghxx ™ can be suppressed. 

4 Conclusions 

In this article we have investigated the viability of regions of large tan j3 in the framework 
of Eat Higgs/A-SUSY models. In the “toy” model we constructed we showed that the 
singlet cubic term is suppressed while the the tadpole and singlino mass parameter terms 
were allowed. Within this framework we showed that there are regions of large tan/3 
that are phenomenologically viable. 

In particular, we computed the one-loop effective potential and showed that the tan [3- 
independent contributions to the Higgs quartic are crucial in raising the Higgs mass to 
the observed value of 125 GeV. We have also shown that non-standard Higgs bosons of 
the same mass as the stops will give comparable contributions to the Higgs quartic when 
A ~ V^Ut- In the region of degenerate non-standard Higgs boson masses the corrections 
are larger than those estimated in Ref. [351 03] . This discrepancy is purely due to the 
assumptions made in Ref. [35l 03] that lead to a split spectrum of heavy CP-even and 
CP-odd scalars. 
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Furthermore, we pointed out that contributions of the neutralinos and charginos to 
electroweak precison observables are small even for large tan /3 when /igs — 500 GeV and 
/i' > 100 GeV. Such large values of /igfj make this region of Fat Higgs/A-SUSY parameter 
space slightly more unnatnral than the low tan /3 region considered in Ref. |22l l23l [2^ 
[25l [26l EZl [28l |29l |30l [SU [32l [33l |3l] . Additionally, this scenario corresponds to the 
decoupling limit where the mixing between the heavy Higgs states and the SM-like 
Higgs is suppressed. Consequently, SM-like Higgs decay properties are with 1% of their 
corresponding Standard Model valnes. Detecting this scenario at the LHC, therefore, 
will be challenging. 

We also found regions of large tan^d in Fat Higgs/A-SUSY models that satisfy all the 
above constraints and provide a viable dark matter candidate. For large n' and tan (5 we 
showed that four possible viable parametric scenarios exist. The dark matter in these 
scenarios could be either most singlino or bino and, depending on their mass, could 
generate the observed relic abnndance through resonant annihilation or co-annihilation. 
In each of these scenarios, direct detection cross-section can be probed at the XENONIT 
experiment. Another possibility is that of a well-tempered neutralino. This scenario 
typically occurs at low values of tan/3 and A^l, where the lightest neutralino’s Higgsino, 
bino and singlino fractions are such that its coupling to the Higgs boson is suppressed. 
The XENONIT experiment may not be able to completely probe this scenario. 
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A Decoupling behavior at one-loop level 


We need to use the tadpoles at the one-loop level to solved for the one-loop corrected 


soft squared mass parameters. Extending Eq. (12) to one-loop order leads to the system 
of three equations. 


Ti = 


dVi 


Higgs 


SI'S. 


d(j)i 


d(t)i 


+ 


5AU 




90* 


= -h ATi = 0, 




i = l,2,3. (42) 


We again can try to solve for the soft masses and in terms of the Higgs 


VEVs. Note that while each as given in Eq. (12), contains only its corresponding 


soft mass ATj in general contain all three soft mass terms. Although obtaining 

the solutions to such a system of equations maybe straightforward, the computation 
conld become complicated when we expand the fnll potential around the trne electroweak 
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symmetry breaking minimum. We can avoid this difficulty by solving Eq. (42) iteratively. 
We first solve for the tree level soft mass squared parameters ’ ("^s) 

using Eq. (12) and then substitute them into ATj. This approximation linearizes Eq. 


which leads to the one-loop corrected soft mass squared parameters solution 


= 


/ 2^0 
= [mi) - 


1 




IGtt^ 

j=D,S 


atil 

dcki 


+ .. 


(43) 


where (m?) is the tree-level solutions of Eq. (12), = bj Vi = {vu,Vd,s), i = 

{H^, S) and (j)i = S). Substituting these solutions into the total potential 

and expanding it about the electroweak symmetry breaking minimum we observe that 
corrections to the CP-even Higgs mass matrix takes the form 


(AM|o)“^ = 


1 


IGtt^ 


E 

i=D,S 


dh\ 


,02 


dcjiadct^h 


1 


db<l 




2U(i d(f)a 


Sab 




Ml, 


+ ..., 


(44) 


By the symmetries of the model, the only field dependences at quadratic order in 6 ? are 


h^, h^, huhd and hj. Thus Eq. (44) suggests that the coefficient of M^^ in the self-energy 


corrections vanishes and that of Mj^^ will be proportional to 


/ 


•5/3 C/3 


V 



(45) 


When these correction are rotated into the basis defined in Eq. (16) we see that the (2, 2) 


element is the only non-zero element. Therefore the decoupling is manifest even at the 
one-loop level. 

B Effective potential derivation 


In this section we apply the procedure outlined in Sec. |2.2| to the computation of one-loop 
radiative corrections from the Higgs sector. First, we deal with degenerate pseudoscalars, 
so that all the one-loop corrections come from a single heavy scale. We will call this Case 
(A). Next, in Case (B), we inspect the effect of splitting the pseudoscalar masses on the 
one-loop corrections, where they now come from two heavy scales. For simplicity, the 
soft terms A\, , m 3 , m^ are taken to vanish in this case. 


(A) Degenerate pseudoscalars 


From the CP-odd mass matrix in Eq. (14), we impose the necessary and sufficient 


condition for mass degeneracy in the pseudoscalars given by ~ 
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(M ^)22 = M\, to obtain 

M = ^A, 

ip = {M^spcp - ml)/\ - 2Axs, 

= -M\s - Ax{M\spcii - ml- \v^S2ii)/\ ( 46 ) 

Respecting this condition, the field-dependent mass matrix for the charged sector is 

“ hi) + ^hl, 

2 

^12 = “ >^^)huhd + 2\Ax{hs - s) M\spcii, 

^22^ = -^{hl- hi) + ^hl, ( 47 ) 

for the CP-odd sector it is 

^11 = + ^^{hl + hi) + ^{hl- hi), 

Mi2 = 2XAx{hs — s) -|- M\sjdCp, 

M 22 = mjj^ + X^{hl + hi) -^{hl- hi), 

Mfa' = 0, 

Mfa' = 0, 

-^33 = ^-1 + -^^(^« + hi) + Af- ml , 

( 48 ) 

and for the CP-even sector it is 

-^n = fTT'Hu + ^'^(hl + hi) + ^{Shl - hi), 

2 

Mi2^ = (2A^ - ^)huhd - 2\Ax{hs - s) + Mls/jC/s, 

M|/ = mjj^ + X'^ihl + hi) - ^{hl - 3hl), 

Mfg^ = 2X^{huhs - Axhd), 

= 2X^{hdhs - Axhu), 

^33 = "^1 + ^^{hl + hi) + A?^ + m'g 

( 49 ) 

The eigenvalues of the charged matrix are given by M ^2 = ^(Ttc T \J^l ~ 4Detc), 
where Trc = and Detc = We only include the 
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_l_2 

contribution from the heavier eigenstate corresponding to ■ Note that when we 
take the supertrace in the charged higgs sector, we obtain a multiplicative factor of 2 
since each charged higgs state comprises of two real physical states. In other words, the 
supertrace is here taken over the full 4x4 squared-mass matrix and not the 2x2 version 
that is usually written down for brevity. 

The eigenvalues of the CP-odd matrix are obtained in a straightforward manner, 
since the upper left 2x2 block is decoupled from . The squared eigenmasses 
are obtained as = 5 (Trp — y^Tr^ — 4Detp), M|p = ^(Tr^ -|- y^Tr^ — 4Detp) and 

M|p = M|3, where Trp = -h and Detp = 

Obtaining the CP-even eigenvalues is non-trivial since we need to deal with a rank 
3 matrix. However, we can take advantage of the degeneracy of the CP-odd scalars by 
employing the following simplifying trick. 

First, consider the characteristic equation of the CP-even matrix, written as 

o o 

asx + a2X + aix -h ag = 0, 


whose solutions are the field-dependent eigenmasses Mf^. The coefficients a,, in terms 
of the matrix elements in Eq. (49), are 


03 

02 

Ol 

oo 


1 , 


-{Mfi + M§2 + ), 

-^U ^22 + ^22 ^33 + ^33 ^11 

o „o „o „o 


-Mi,\ 






O 2 


We also know, in terms of the eigenmasses, that 

02 = —{Mfg + -|- M^g), 

Ol = M^Ml, + + Ml,Ml 


(51) 


Now the CP-even sector contribution to the effective potential, from Eq. (18), is 

1 


The quantity in brackets can be re-written using Eq. © as simply 


+ (Ml^r = «2 - 2 ai - 


(52) 


(53) 


The coefficients oi and 02 may be read off Eq. (50), while we may still have to determine 
Mf g analytically. This is, however, a simple task if we write ^ as a power 
Ml: 


senes in 


= 4 + 1 ), 


(54) 
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where bi is at most quadratic in the background fields. Putting Eqs. and 54 into 
Eq. (52), we obtain the one-loop effective potential contribution simply as 


^ ^- 2ai - bj] log , (55) 

where we have discarded 0(1/M^) terms that are irrelevant in obtaining the required 
self-energy corrections. 

After including all the one-loop corrections, the final expressions we obtain for the 
CP-even mass matrix are now as follows. 


+ Hii; 

22 ^ + ^ 22 ) 
23 ~ ^23; 

where 


(m^) = (2A\2 - m| - mI)s^C0 + U 12 ; 

= 2Au^eff'S/3 + His; 

(M^)^^ = M^ + n33, 



Ml 



log 



(56) 





and 


nil = ^^[-32AS2(2 c2;3 - sip) + 2AV(3c 2/3 - l)(3si^ + 2) 
+5'^(4c^ + — 7 sip — 1 — C2/3(4c^ — 4c^ + 5^2^ + 3)) 

+64A2^^^cot/3]log(^^^ , 

ni2 = 25071^ ~ 2) — 2A^g'^S2/3(8c^ — ISs^^ + 14) 

+5^'S2/3(4c^ + 4c^ — 7 sip + 3) 


-64A' 




[32A^c^(2c2;3 + sip) - 2AV(3 c2P + l)(3si^ + 2) 

„4 ,1 „2 I '7„2 I 1 I „ _/'/i„4 j 2 I c „2 I 


256712 

+5^(4c)^ — 4c^ + 7^2^ + 1 + C2/3(4 c^ — 4c^ + + 3)) 

+64A2^^^ tan/3] log , 


His 


n- 


-23 


= ^*^[ 12 A^s^ + Xg^sp{3c2fs + 2 c^ + 1 ) 

/^eff 

+ ^^ 32772^^ ^ 

= ^^[ 12 A^c^ + Xg‘^cp{-3c2p + 2 c^ + 1 )) 

Meff 

+ ~' ^ 7 r™^ ' ^ ^ ^ + ^C 2 P + 3cip)'^]log , 

/M2 a 

+A(64s^ + 29v‘^)s20 + Xv'^sep) + g‘^v^S2g{3 + 4c^ + C 4 ^)]}log j i^'^) 

If we set all NMSSM-specific parameters to zero in the above, we recover the MSSM 
limit presented in [601 in EH 163] , The soft term A\ decouples at one-loop order and 
does not contribute to the SM Higgs quartic coupling, a property best seen in the basis 
of Eq. (16). The SM Higgs boson mass is then identified as 

^hh = ^^v^slp + M^clp -b Hhh, 


V 


H/i/i — g227r2 ^ 04/3 — 3csp) -\- 4A^(7^(—9 — 4c^ + (4c^ — 2 ) 04/3 + ^cs/s) 

—g'^(—11 + 8c^ — 16 c^ -|- 8c^C4/3 -|- 3c 8/3)] log ] ( 58 ) 
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A\ is absent in the expression above, confirming its decoupling behavior at the one- 
loop level. Moreover, if we neglect the electroweak strength corrections, in the limit of 
large tan /3 we get 


TT 

hm Uhh = -r^ 
tan/3>l 


, . Ml 


in agreement with our heuristic estimate in Eq. (25). 


(59) 


(B) Non-degenerate pseudoscalars: a simple case. 


We now show the effect of a split pseudoscalar spectrum on the radiative corrections. 
For simplicity, we assume the parameters A\, fi',ms,m!g vanish. We also neglect g- 
dependent terms in the one-loop piece, since the largest contributions to the SM Higgs 
quartic in our model arise from the A-dependent terms. With these simplifications, the 
field-dependent squared mass matrices for the charged, CP-odd and CP-even sectors are 
respectively given by 


+ A^h^, Mf/ = X^Khd -h Ml^spCfi, 


= + (60) 


M,f = m%^ + X\hl + hl), 

uff = Ml^sfjcp, 

M^2=ml^ + X\hl + hl), 

Mfs" = 0, 

Mfg' = 0, 

M[,^ = ml + X\hl + hl) 


(61) 


+ A^(h^ + hg), 

AI 22 = + A^(/i^-I-/ig), 

AI 2 Z — ‘^X'^hdhs, 


M^2 = ‘^>^'^huhd - M^^sgcid, 
= 2X^huhs, 

Mi,^ = ml + X\hl + hly, 


(62) 


Obtaining the eigenvalues of the charged and CP-odd systems is straightforward 
again, as we found in Case (A). To obtain the eigenvalues of the CP-even matrix, we 
solve for the roots of its characteristic equation (a cubic polynomial) as a power series 
in and 

After collecting the one-loop contributions from all three sectors and summing over 
them, we obtain the CP-even mass matrix as 

(m^)+ Hn; (m^) = {2X\^ - M| - + His; 



M|c^ + + II22; 

= 2XvgeSCl3 + n23; 
23 



‘^Xvfiessg + Hia; 



— + Has 
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where 


MAn = 



/^eff 


8 vr 2 Ml^ - Ml 


log 


'Ml 


As 




Mas 


= M 




(63) 


and 


nil = 


ni 2 = 


AS 2 2 
IGvr^ ^ 


-{4c., + c., + l)log(^]+21„gg 


^W (2 + C4,)log|:^ 


Ml 


4., 2 


n ^ ^ ^ 
A^n/ieff 


(Ml \ [Ml' 

-(-4c2/3 + C 4;3 + 1) log ( I + 2 log ' ® 


His = 


n23 = 


Has = 


IGvr^ 

A^n//eff 

IGvr^ 

IGvr^ 


^/3 


M| 


-(l + 3.,,)l„g(^)+41og(^ 


M| 


-Cp 


-(1 - 3C2/3) log 


Ajl 

Ml 


+ 41og 


'Ml 

M| 


log 


Ml 


(64) 

We make the following observations concerning the above expressions. First, notice 
that in the limit Maj^ = Ma^, they are consistent with the results in Case (A) with 
g, Ax -A 0. Second, we observe that corrections from the heavy doublet Higgses are 
/3-dependent and those from the heavy singlet Higgses are not, as reflected in the co¬ 
efficients of log(M^^/M^) and log(M^^/M^) respectively. Third, there is a marked 
difference in contributions from the scales Ma^^ and Ma^ to the SM Higgs quartic, 
which can be understood in the basis of Eq. (16). Rotating Yiij into this basis, the 


correction to the SM Higgs boson mass is identified as 


H/i/i = 


IGtt^ 


(c^(2 -(- C 4 ^) — S^(l C 4 ^ -|- 4c 2/3)) log 


All, 

Ml 


+ 2s^ log 


'yh' 

Ml 


The difference in the co-efficients of the logarithms are greatest at tan /3 
at tan /3 S> 1. In the latter limit, we obtain 


1, and smallest 


TT A^n^ 

hm Yihh = 
tan/3>l IGtT"^ 


2 log 


'Ml \ [Ml' 

40 I +21og I 4o 


Mi 


Mi 


which is consistent with our qualitative estimate in Eq. (24). 
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